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The ground state and zero-temperature magnetization process of the spin-1/2 Ising-Heisenberg 
model on two-dimensional triangles-in-triangles lattices is exactly calculated using eigenstates of the 
smallest commuting spin clusters. Our ground-state analysis of the investigated classical-quantum 
spin model reveals three unconventional dimerized or trimerized quantum ground states besides two 
classical ground states. It is demonstrated that the spin frustration is responsible for a variety of 
magnetization scenarios with up to three or four intermediate magnetization plateaus of either quan- 
tum or classical nature. The exact analytical results for the Ising-Heisenberg model are confronted 
with the corresponding results for the purely quantum Heisenberg model, which were obtained by 
numerical exact diagonalizations based on the Lanczos algorithm for finite-size spin clusters of 24 
and 21 sites, respectively. It is shown that the zero-temperature magnetization process of both 
models is quite reminiscent and hence, one may obtain some insight into the ground states of the 
quantum Heisenberg model from the rigorous results for the Ising-Heisenberg model even though 
exact ground states for the Ising-Heisenberg model do not represent true ground states for the pure 
quantum Heisenberg model. 
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I. INTRODUCTION 

Two-dimensional frustrated quantum spin systems 
have attracted a great amount of research interest over 
the last few decades X Frustration as a phenomenon orig- 
inating from competing interactions incompatible with a 
lattice geometry prevents spins from simultaneously sat- 
isfying all pair spin-spin interactions and may thus lead 
to a variety of fascinating quantum phenomena. Among 
the most interesting properties emerging in frustrated 
quantum spin systems one could mention the existence of 
exotic ground states with a spin-liquid or valence-bond- 
solid character^ - — At sufficiently low temperatures, the 
external magnetic field may additionally cause the ap- 
pearance of intermediate magnetization plateaus at some 
fractional values of the saturation magnetization^— 

SrCu2 (1303)2 represents one of the most famous ex- 
perimental realizations of the frustrated quantum spin 
systems, which display an intriguing sequence of mag- 
netization plateaus before reaching the saturation mag- 
netization. High-field measurements performed on this 
layered magnet have revealed intermediate plateaus at 
|, |, and i of the saturation magnetization^ - — Sub- 
sequently, torque measurements by Sebastian et alM- 
have suggested the presence of additional striking frac- 
tional plateaus. A similar conclusion has also been 
reached more recently on the basis of torque and NMR 
measurements^, with a different sequence of plateaus 
however. From the theoretical point of view, deep in- 
sight into the low-temperature magnetization process of 
SrCu2 (B03)2 has been provided by the investigation of 
the spin-i quantum Heisenberg model on the Shastry- 



Sutherland lattice^—, even if the actual sequence of 
plateaus for large inter-dimer coupling is still debated. 

The existence of an intermediate magnetization 
plateau at one third of the saturation magnetization has 
been recently reported for the family of polymeric coordi- 
nation compounds Cu9X 2 (cpa)6-nH 2 (X=F,Cl,Br and 
cpa=carboxypentonic acid). 2 ? -25 It provides a remark- 
able experimental realization of the spin-i Heisenberg 
model on a quite exotic triangulated kagome lattice^ - — 
It actually turns out that the magnetic structure of this 
series of magnetic compounds is formed by two inequiva- 
lent lattice positions of the spin-i Cu 2+ ions, whereas 
smaller equilateral triangles of Cu 2+ ions (Heisenberg 
trimers) are mutually inter-connected through Cu 2+ ions 
(Heisenberg monomers) situated at lattice points of a 
simple kagome lattice. From the geometric point of view, 
the overall magnetic structure thus consists of smaller tri- 
angles of magnetic ions embedded in the larger triangles 
of a kagome lattice and hence, it belongs to a class of two- 
dimensional triangles-in-triangles (TIT) lattices. Recent 
exact solution for the spin-i Ising-Heisenberg model on 
two closely related TIT lattices descending from a sim- 
ple triangular lattice rather than from a kagome lattice 
gave evidence of a significant impact of local quantum 
fluctuations on the overall magnetic behavior of this hy- 
brid classical-quantum spin models The main purpose 
of this work is to find possible ground states of the spin- 
■i Ising-Heisenberg model in the presence of an external 
magnetic field and to compare its zero-temperature mag- 
netization process with that of the full quantum version 
of the spin-^ Heisenberg model. 

The outline of this paper is as follows. In Sec. |H]we m_ 
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troduce the spin-i Ising-Heisenberg and Heisenberg mod- 
els on two TIT lattices and briefly describe the procedure 
used for the calculation of the zero-temperature mag- 
netization process. The most interesting results for the 
ground state and zero-temperature magnetization curves 
are presented in Sec. HID Finally, the summary of the 
most important scientific achievements and future out- 
looks are mentioned in Sec. IIV1 

II. MODEL AND ITS GROUND STATE 

Let us introduce first the underlying magnetic lat- 
tice for frustrated quantum spin models, whose zero- 
temperature magnetization process will be explored in 
detail. Fig. [T] schematically illustrates two particular ex- 
amples of the TIT lattices, which can be derived from a 
simple triangular lattice by placing an additional inner 
triangle either into each up-pointing triangle (Fig. QJi) or 
into each triangle (Fig. [Tb) of a triangular lattice. The 
two displayed TIT lattices can alternatively be viewed 
as being composed of inter-connected stars (see Fig. [5]), 
which have an equilateral triangle in their center and are 
joined together through outer isosceles triangles attached 
to each side of the inner equilateral triangles. Clearly, 
there are two inequivalent lattice sites in the TIT lattices, 
the ones forming the inner equilateral triangles shown by 
empty circles and the other ones solely belonging to the 
outer isosceles triangles depicted by filled circles. The 
more general anisotropic version of the spin- ^ Heisenberg 
model on two different but topologically related TIT lat- 
tices can be defined through the Hamiltonian 

o—o o—m 

u = j t j2Si-Sj+j s j2 (s k -s l ) A -hjy%, (i) 

<i,j> <k,l> p 

where (S fe • S ; )a = A(5^Sf + S y k Sf) + S z k S[ and S< = 
(Sf , Sf , Sf) denotes the standard spin-i operator. The 
interaction term J t describes the isotropic Heisenberg in- 
teraction between the nearest-neighbor spins from the in- 
ner equilateral triangles while the interaction term J s la- 
bels the anisotropic XXZ Heisenberg interaction between 
the nearest-neighbor spins from the inner equilateral and 
outer isosceles triangles, respectively. Finally, the last 
Zeeman's term h accounts for the magnetostatic energy 
of all spins in an external magnetic field. In what follows, 
our main attention will be concentrated on two special 
cases of the Hamiltonian ((T|) either with A = or 1, 
respectively. While the Hamiltonian ([1]) with A = re- 
duces to the hybrid classical-quantum Ising-Heisenberg 
model^ the Hamiltonian ([1]) with A = 1 corresponds to 
the isotropic quantum Heisenberg model. 

A. Ising-Heisenberg model 

First, we turn our attention to the spin-^ Ising- 
Heisenberg model defined through the Hamiltonian ([1]) 




FIG. 1: Two particular examples of the TIT lattices derived 
from a simple triangular lattice. The lattice positions solely 
belonging to the inner equilateral (outer isosceles) triangles 
are schematically represented by empty (filled) circles. 




FIG. 2: Elementary spin cluster in the form of a spin star. 

with A = 0. It is worthwhile to remark that the spe- 
cific choice A = disregards the transverse part of 
the exchange interaction J s , which consequently becomes 
the pure Ising coupling. Owing to this fact, the spins 
from the outer isosceles triangles are merely coupled to 
their nearest neighbors by the Ising interaction, while the 
isotropic Heisenberg interaction J t between the nearest- 
neighbor spins from the inner equilateral triangles is still 
preserved. For easy reference, the former spins from the 
outer isosceles triangles will be therefore referred to as 
the Ising spins and the latter spins from the inner equi- 
lateral triangles as the Heisenberg spins. 

For further convenience, the total Hamiltonian ([1} can 
be decomposed into a sum of cluster Hamiltonians % a 

-yN 

U = j2Ua, (2) 

a=l 

where N is the total number of Ising spins, 7 = 1 (7 = 2) 
for the first (second) TIT lattice, so that 7^ is the total 
number of six-spin clusters schematically visualized in 
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Fig. [2] The cluster Hamiltonian 7i a is given by 



D J 

'fL a = Jt ^<M ■ S a ,j+1 + J s ^( Sg,i^a,i+3 + Sa,i&a, 



i+4 



i=4 
6 



»* — /I ' ^ — 1 



with the convention 5a f 7 = S^. 

It is worth noticing that the factor 1 



(3) 



the 



3 7 emerging in 
last term ensures correct counting of the Zeeman's en- 
ergy of the Ising spins, which is equally split into the 
37 different cluster Hamiltonians involving a given Ising 
spin. Moreover, it is of fundamental importance that 
different cluster Hamiltonians commute with each other 
[Hi^T-Lj] = as they have in common at most one Ising 
spin. Accordingly, different cluster Hamiltonians can be 
diagonalized independently of each other as they belong 
to mutually orthogonal Hilbert subspaces and one con- 
sequently gets the factorized ground state as a tensor 
product over the lowest-energy eigenvectors of the clus- 
ter Hamiltonians ([3]). The minimization of all clusters 
simultaneously is not always possible if they share com- 
mon spins, but in our case the different ground-state ver- 
ify this condition. 

In order to hnd the exact ground state of the spin-^ 
Ising-Heisenberg model on the two considered TIT lat- 
tices, it is thus sufficient to solve the eigenvalue problem 
for the cluster Hamiltonian ^ by considering all possi- 
ble configurations of three Ising spins involved therein. 
For the sake of brevity, let us merely list all eigenval- 
ues of the cluster Hamiltonian ([3]) without quoting the 
corresponding eigenvectors, which will be later specified 
just for the eigenstates that may become the ground 
state. The eigenvectors of the cluster Hamiltonian can 
be classified by the total magnetization of the Ising spins, 
because this operator commutes with the Hamitonian. 
The eigenvalues of the cluster Hamiltonian ([3]) for three 
Ising spins equally aligned into the external-field direc- 
tion ISMoSfa) = |ttt) read 
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The energy spectrum of the cluster Hamiltonian (j3|) for 
three Ising spins aligned opposite to the external mag- 
netic field (\S* 2^0,3) = I III)) can be obtained by 
changing the sign of the magnetic field h in Eq. (|4]). 

If two Ising spins have the same orientation with re- 
spect to the external magnetic field and one Ising spins 
is in the opposite direction, for instance \S* iS* 2$! 3) = 



ITTI) 1 one obtains the following eigenvalues of the cluster 
Hamiltonian 
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We can get the energy spectrum of the cluster Hamil- 
tonian ([3J for one Ising spin pointing in the direction of 
the external magnetic field and two Ising spins aligned in 
the opposite direction \S* ^Sa^Sa^) = Itll) by changing 
the sign of the magnetic field h in Eq. 

The lowest-energy eigenvalue from the set given by 
Eqs. (j4])-([5|) unambiguously determines the ground state 
of the Ising-Heisenberg model, which will be investigated 
in detail in Sec. |Hl] depending on the relative strength of 
both interaction constants and the magnetic field. In ad- 
dition, the zero-temperature magnetization can easily be 
obtained from the ground-state energy using the relation 



27 dE g 
'37+I dh 



(7) 



where the factor 



27 
37+1 



ensures the proper normalization 



of the total magnetization with respect to its saturation 
value since the elementary unit cell consists of one Ising 
spin and three (six) Heisenberg spins for the first (second) 
TIT lattice displayed in Fig. [TJi (Fig. [Tb). 



B. Heisenberg model 

The other particular case of the Hamiltonian ([1]) with 
A = I corresponds to the isotropic spin-i Heisenberg 
model in a presence of the external magnetic field. Of 
course, the ground state of the fully quantum Heisenberg 
model on the TIT lattices cannot be rigorously obtained 
by analytical calculations and one has to rely on numeri- 
cal methods. The lowest-energy eigenstates of the spin-i 
Heisenberg model on two TIT lattices were obtained by 
employing the numerical exact diagonalization based on 
the Lanczos algorithm. Our numerical calculations were 
limited to relatively small finite-size subsystems when 
comparing them with the relevant magnetic unit cell of 
the classical ground state of both TIT lattices, which ex- 
tends over three lattice unit cells. The magnetic unit cell 
of the classical ground state for the first TIT lattice thus 
contains 12 spins (the elementary unit cell consists of 
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FIG. 3: (a) The magnetic unit cell of the first TIT lattice is 
shown as the 12-site rhombus spin cluster delimited by thick 
(orange) line, which also includes spin orientations within the 
classical ground state. The 24-site parallelogram spin cluster 
constituted by two magnetic unit cells was used for the exact 
diagonalization as the largest spin cluster; (b) The magnetic 
unit cell for the second TIT lattice is formed by the 21-site 
rhombus spin cluster, which was used for the exact diagonal- 
ization as the largest spin cluster (spin orientations within the 
classical ground state are not shown for clarity). 



4 spins), while our exact numerical diagonalization was 
limited to the rhombus and parallelogram spin clusters 
of 12 and 24 sites as specified in Fig. On the other 
hand, the magnetic unit cell of the classical ground state 
for the second TIT lattice already contains 21 spins (the 
elementary unit cell consists of 7 spins), which was also 
an upper limit for the spin-cluster size of our numerical 
diagonalization (see Fig. [3Jd). 

Obviously, the total spin operator S* ot — J2 P com- 
mutes with the Hamiltonian [H, S^ ot ] = 0, which means 
that the external field does not alter the relevant eigen- 
states and the total magnetization S* at = J2 P Sp is a 
conserved quantity. Hence, the full energy spectrum of 
the spin-| Heisenberg model in the presence of the ex- 
ternal magnetic field can be obtained from the energy 
spectrum in a zero magnetic field according to 



E i (h,S? ot ) = E i (0,S? ot )-hSZ 



(8) 



It is quite clear from Eq. © that the ground state may 
correspond only to the lowest-energy eigenstate from 
some sector with the fixed value of Sf ot . The field- 
induced change in the ground-state spin arrangement 
can thus be related to crossings between the lowest- 
energy levels from the sectors with different S^ ot . Hence, 
the zero-temperature magnetization curves can simply 



be constructed from the lowest-energy eigenstates, which 
were obtained with the help of numerical exact diagonal- 
ization based on the Lanczos algorithm for all possible 
sectors of Sf ot at zero magnetic field. 



III. RESULTS AND DISCUSSION 

Before proceeding to the discussion of the most in- 
teresting results for the zero-temperature magnetization 
curves of the spin-i Ising- Heisenberg and Heisenberg 
model on the TIT lattices, let us introduce a more con- 
venient parametrization of the interaction constants 



Jt = J sin a, J s = J cos a, 



a e (0,90 c 



(9) 



which is useful in an attempt to provide a comprehensive 
analysis of the full parameter space. The specific choice 
a = corresponds to Jt = and J s = J, which means 
that the isotropic Heisenberg interaction along the inner 
equilateral triangles vanishes and the nearest-neighbor 
spins from the inner equilateral and outer isosceles trian- 
gles are merely coupled through the interaction J s = J. 
On the other hand, the opposite limit a = 90° im- 
plies that the investigated spin systems break into a set 
of independent Heisenberg trimers and monomers since 
Jt = J and J s — 0. For simplicity, our subsequent anal- 
ysis will be restricted just to the most interesting partic- 
ular case with the antiferromagnetic interactions Jt > 
and J s > 0, whose relative strength is determined by the 
newly introduced parameter a through Jt/J s = tana. 
Last but not least, it should be pointed out that the an- 
tiferromagnetic Heisenberg interaction J t is responsible 
for a spin frustration, which becomes relevant for suffi- 
ciently high values of the interaction ratio Jt/J s (or a). 

For the sake of comparison, the zero-temperature mag- 
netizations of the spin-i Ising- Heisenberg and Heisenberg 
model on the first TIT lattice are plotted in Fig. |4] against 
the magnetic field H/J and the parameter a. First, let 
us make a few comments on the magnetization curves 
of the spin-i Ising- Heisenberg model for which the re- 
spective ground states can be found quite rigorously. As 
one can see from Fig. H^,, the Ising- Heisenberg model dis- 
plays five possible ground states two of which are classical 
and three quantum in character. If the interaction ratio 
Jt/Js < 3 (i.e. a < 33.7°), the classical ferrimagnetic 
state emerges in the ground state at low enough fields as 
a result of the following lowest-energy eigenstate of the 
spin star (see Fig. [2] for the definition of the sites inside 
a star) 



1 1) = \lll)al.a2,a3® | ttt)o4,o5,o6, 

3 r 3 r 3 , 1 . 

-E/t = —Jt j s n H h, 

4 2 2 2 7 ' 



(10) 



which indicates that the Heisenberg spins are fully 
aligned towards the magnetic field in contrast to the 
Ising spins pointing in an opposite direction (see Fig. [S^ 
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FIG. 4: (a) A three-dimensional plot for the zero-temperature 
magnetization of the spin-i Ising-Heisenberg model on the 
first TIT lattice (Fig. \T^,) as a function of the magnetic field 
H/J and the parameter a; (b) The same three-dimensional 
plot but for the pure quantum spin-i Heisenberg model on 
the 24-site parallelogram spin cluster depicted in Fig. 



for a schematic illustration of this classical spin arrange- 
ment). The total magnetization of the classical ferrimag- 
netic state will be consequently equal to one half of the 
saturation magnetization, because the magnetization of 
one Ising spin just partially compensates the magnetiza- 
tion of three Heisenberg spins per elementary unit cell 
of the first TIT lattice. The classical ferrimagnetic state 
ends up at the saturation field h s — 3J S , above which the 
saturated paramagnetic state with a perfect alignment of 
all spins into the external magnetic field emerges owing 
to the lowest-energy eigenstate (see Fig. [5p) 



|H) = |ttt)al,a2,a3® | Ttt)a4,a5,a6, 



tt-^s — -^h — —h. 
2 2 2 7 



(11) 



formation of an intermediate magnetization plateau at 
half of the saturation magnetization comes from the dom- 
inant character of the Ising interaction J s that governs 
the magnetic behavior in the relevant parameter region. 

A more striking situation appears for the other pos- 
sible choice of the interaction ratio Jt/J s > § (i-c 
a > 33.7°), which promotes spin frustration and the 
gradual strengthening of local quantum fluctuations due 
to the crucial role of the Heisenberg interaction J t . Under 
this condition, the notable up-down-down (udd) dimer- 
ized ground state can be detected at low enough fields 
that can be constructed from the eigenvector 
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(12) 



The spin alignment inherent to the udd dimerized state 
unambiguously given by the eigenstate (|12[) is schemat- 
ically depicted in Fig. [5};. It can be easily understood 
from this figure that one out of the three Ising spins be- 
longing to the same spin star (Fig. [5]) is aligned towards 
the magnetic field, while the other two are oriented in 
an opposite direction. In addition, there is exactly one 
field-aligned Heisenberg spin and one singlet dimer per 
equilateral triangle of the Heisenberg spins (Heisenberg 
trimer), whereas up-pointing Ising spins are surrounded 
by singlet dimers and up-pointing Heisenberg spins by 
down-pointing Ising spins. Altogether, the total mag- 
netization i per Heisenberg trimer is just partially com- 
pensated by the average magnetization — ^ per Ising spin 
and hence, the udd dimerized state (|12[) manifests it- 
self through the magnetization plateau at one sixth of 
the saturation magnetization. If Jt/J s & (I'fv^— §) 
[or equivalently a £ (33.7°, 43.6°)], the magnetization 
curves show two intermediate magnetization plateaus at 
one sixth and one half of the saturation magnetization, 
which result from two consecutive field-induced transi- 
tions between the udd dimerized state, the classical fer- 
rimagnetic state and the saturated paramagnetic state. 

Besides, one also encounters an additional magnetiza- 
tion plateau at one third of the saturation magnetiza- 
tion for stronger values of the interaction ratio J t /J s > 
I V7 — |, which reflects an appearance of another intrigu- 
ing up-up-down (uud) trimerized ground state stemming 
from the lowest-energy eigenstate 
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As could be expected, this classical mechanism for the 



where the mixing angle (f> for a quantum superposition 
of three uud states of the Heisenberg trimers is given by 

. A fragment from the spin config- 



= ^ arctan 



+2J S 



6 

<ZE§> = (|w> - l-t-t»/V2 




(a) classical ferrimagnetic state (I) (b) saturated paramagnetic state (II) (c) udd dimerized state (III) 




(d) uud trimerized state (IV) (e) uuu dimerized state (V) 



FIG. 5: Two classical and three quantum ground states of the spin-| Ising-Heisenberg model referred to as: (a) the classical 
ferrimagnetic state I; (b) the saturated paramagnetic state II; (c) the uud dimerized state III; (d) the udd trimerized state IV; 
(e) the uuu dimerized state V. 



uration of the uud trimerized state is schematically illus- 
trated in Fig. [5jl. Apparently, there are two up-pointing 
and one down-pointing Ising spin per elementary spin 
star (Fig. [5]) and this non-uniform spin alignment of the 
Ising spins is also transfered to a quantum superposi- 
tion of three uud states of the Heisenberg trimers when 
the highest occurrence probability for a down-pointing 
spin belongs to the Heisenberg spin attached to the two 
up-pointing Ising spins. The total magnetization of the 
uud trimerized state will be determined by the magne- 
tizations i per Ising spin and \ per Heisenberg trimer, 
which will finally lead to the other magnetization plateau 
at one third of the saturation magnetization. Bearing 
this in mind, the magnetization curves exhibit three in- 
termediate magnetization plateaus at one sixth, one third 
and one half of the saturation magnetization, which are 
consistent with three consecutive field-induced transi- 
tions between the udd dimerized state, the uud trimer- 
ized state, the classical ferrimagnetic state and the satu- 
rated paramagnetic state whenever Jt/J s & (f V7— §, §) 
[a e (43.6°, 53.1°)]. 

If the interaction ratio exceeds Jt/J s > |, there ap- 
pears another unusual quantum up-up-up (uuu) dimer- 
ized ground state arising out from the following eigen- 
states 
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It is noteworthy that all the Ising spins are fully aligned 
in the uuu dimerized state towards the magnetic field, 
whereas the Heisenberg trimers consist of one singlet 
dimer and one field-aligned Heisenberg spin as depicted 
Fig. [5^. Hence, it follows that the uuu dimerized state 
is highly macroscopically degenerate in contrast to all 
the previously described ground states, because two lin- 
early independent eigenstates can be constructed for each 
Heisenberg trimer from the inner equilateral triangles. 
Even though the magnetization curves for J t / J s > | will 
still contain three intermediate plateaus at one sixth, one 
third and one half of the saturation magnetization, the 
last magnetization plateau before the saturation now cor- 
responds to the uuu dimerized state rather than to the 
classical ferrimagnetic state. 

Next, let us compare the magnetization process of 
the Ising-Heisenberg model (Fig. |4ji) with the relevant 
numerical data for the analogous but purely quantum 
Heisenberg model (Fig. |4Jd). Except for a few mi- 
nor differences to be specified in the following, the 
zero-temperature magnetization curves of the Ising- 
Heisenberg model and of the quantum Heisenberg model 
have several important common features. In particu- 
lar, the plateaus identified at one sixth, one third and 
one half of the saturation magnetization in the Ising- 
Heisenberg model can be clearly recognized in the mag- 
netization curve of the Heisenberg model as the largest 
steps, the smaller ones being very likely finite-size ef- 
fects for regions in which the magnetization increases 
smoothly for the infinite system. Of course, it would be 
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necessary to have several sizes at hand and to perform 
a finite-size analysis to be completely sure about the ac- 
tual sequence of the plateaus in the Heisenberg model, 
but the remarkable agreement with the Ising-Heisenberg 
model makes it likely that all the plateaus of the Ising- 
Heisenberg case survive. Thus, the only significant differ- 
ence in the magnetization process of the Ising-Heisenberg 
and Heisenberg models appears to be the presence of di- 
rect magnetization jumps between plateaus in the former 
classical-quantum Ising-Heisenberg model, while gapless 
phases in which the magnetization increases smoothly 
between plateaus are expected to be generically present 
in the Heisenberg model. For instance, it is quite clear 
from Fig. |4Jd that the magnetization plateau at one half 
of the saturation magnetization in the parameter region 
J s 3> Jt is followed by a smooth continuous increase of 
the magnetization over a relatively wide field region be- 
fore reaching the saturation magnetization. 

Finally, let us bring some insight into the nature of the 
gapful phases of the quantum Heisenberg model, which 
correspond to the most robust magnetization plateaus 
observable at one sixth, one third and one half of the 
saturation magnetization. Even though the exact ground 
states of the Ising-Heisenberg model do not represent the 
true ground states for the purely quantum Heisenberg 
model, they might provide at least some useful hint to 
the true ground states of the quantum Heisenberg model. 
It should be stressed, however, that the spins from the 
outer isosceles triangles do not represent an insurmount- 
able barrier for quantum fluctuations in the Heisenberg 
model, which may consequently cause a quantum reduc- 
tion of the magnetization for both kinds of spins from 
the inner equilateral as well as outer isosceles triangles. 

The nature of the ground states is most evident for 
the two gapped phases which lead to the presence of 
an intermediate one-half plateau in the magnetization 
process. The first one-half plateau that emerges in the 
region with dominant interaction J s corresponds to a 
quantum ferrimagnetic phase, a quantum analog of the 
classical ferrimagnetic state shown in Fig. [5^,. Our ex- 
act diagonalization data for the nearest-neighbor spin- 
spin correlation functions of the largest 24-site parallelo- 
gram spin cluster indeed provide a strong support of this 
statement. Considering for instance the particular case 
with Jt/J s = (i.e. a — 0°), the correlation functions 
(S 4i _ 2 'S/ii-i) =0.23256 and (S 4i _i ■ S 4l ) = (S 4i • S 4l _ 2 ) = 
0.23385 imply a strong ferromagnetic correlation between 
the nearest-neighbor spins from the inner equilateral tri- 
angles, see Fig. [3^ for the relevant numbering of lattice 
sites to be congruent modulo 24. By contrast, the corre- 
lation functions (S 4 j_3 • S 4 j) = (S 4 ^_3 • S 4 i + i6) = —0.32053 
and (S 4 i_3 • S 4 i_5) = (S 4 ^_3 • S 4 j_g) = (S 4 i_3 • S 4 ;_2) = 
(S 4i _ 3 • S 4i _ 9 ) = —0.32825 indicate a relatively strong an- 
tifcrromagnctic correlation between the nearest-neighbor 
spins from the inner equilateral and outer isosceles trian- 
gles, respectively. On the basis of these results, it could 
be concluded that the one-half magnetization plateau in 
the region with dominant interaction J s indeed corre- 



sponds to the quantum ferrimagnetic phase in which both 
kinds of magnetic moments are subject to a quantum re- 
duction of the magnetization in contrast to the fully sat- 
urated magnetic moments of the classical ferrimagnetic 
state of the Ising-Heisenberg model. Note that the small 
differences in the relative strength of the nearest-neighbor 
pair correlation functions can be attributed to the asym- 
metry of the 24-site parallelogram spin cluster. If the 
Heisenberg interaction J t inside the inner equilateral tri- 
angles is turned on, the quantum reduction of the mag- 
netization is enhanced due to the spin frustration origi- 
nating from this interaction term. As a matter of fact, 
one generally observes a gradual weakening of the corre- 
lation functions between the nearest-neighbor spins from 
the inner equilateral triangles (S 4 j_2 ■ S 4 i_i) =0.22449, 
(S 4l _i • Sm) = (S 4i • S 4i _ 2 ) = 0.22671, as well as, the 
correlation functions between the nearest-neighbor spins 
from the inner equilateral and outer isosceles triangles 
(S 4 j_3 • S 4 j) = (S 4i _ 3 • S 4i+ i 6 ) = —0.31836, (S 4 i_3 • 
S 4 i— 5) = \S4i-3 • S 4 i_6) = —0.32675, (S 4 j_3 • S 4 i_2) = 
(S 4 i_3 • S 4 ;_g) = —0.32886 as evidenced by the above 
results calculated for the particular case with Jt/J s = 5 
(i.e. a = 26.6°). 

Next, let us discuss the nature of the other one-half 
magnetization plateau, which appears in the region with 
dominant interaction J t . Let us recall first that the ex- 
act ground state of the Ising-Heisenberg model in this 
parameter range is constituted by the macroscopically 
degenerate uuu dimerized state. The uuu dimerized 
state includes the fully polarized Ising spins from the 
outer isosceles triangles, while there is one singlet dimcr 
and one polarized spin per Heisenberg trimer residing 
on the inner equilateral triangle (see Fig. [SJs). Due to 
the macroscopic degeneracy of the uuu dimerized state, 
one should also expect a highly resonating character of 
the relevant ground-state manifold with the singlet-dimer 
state being equally distributed over all three bonds of 
each inner equilateral triangle. Based on these consid- 
erations, the correlation function between the nearest- 
neighbor Heisenberg spins from the inner equilateral tri- 
angles should be equal to (S, • Sj) = —0.25, whereas 
the correlation function between the nearest-neighbor 
Heisenberg and Ising spins from the inner equilateral and 
outer isosceles triangles should be equal to (S%Sf) = ~ 
0.08333. The exact diagonalization data for the nearest- 
neighbor pair correlation functions bear evidence that the 
ground state of the pure quantum Heisenberg model is 
quite reminiscent of the uuu dimerized ground state of 
the Ising-Heisenberg model. For illustration, let us quote 
the correlation functions between the nearest-neighbor 
spins from the inner equilateral triangles of the 24-site 
parallelogram spin cluster (S 4i _ 2 • S 4i _i) = —0.24500, 
(S 4l _i • S 4i ) = (S 4i ■ S 4l _ 2 ) = -0.24431 along with the 
correlation functions between the nearest-neighbor spins 
from the inner equilateral and outer isosceles triangles 
(S 4 i_3 • S 4 j) = (S 4 i_3 • S 4 i + i6) =0.06370, (S 4 i_3 • S 4 i_5) = 
(S 4 i_3 • S 4 i_6) =0.06366, (S 4i _ 3 • S 4i _ 2 ) = (S 4i _ 3 ■ S 4i _ g ) = 
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0.06397 calculated for the particular case with Jt/J s = 4 
(i.e. a — 76°). It could be easily understood from 
these results that the true ground state of the Heisen- 
berg model quite closely resembles the uuu dimerized 
state of the Ising-Heisenberg model, which means that 
the magnetic moments of the Heisenberg spins from the 
outer isosceles triangles underlie only a small quantum re- 
duction of the magnetization (they are almost fully sat- 
urated) in contrast with the magnetic moments of the 
Heisenberg spins from the inner equilateral triangles un- 
derlying a quantum reduction of the magnetization ap- 
proximately to one third of the saturation value. As far 
as the other two gapped phases corresponding to the one- 
sixth and one-third plateaus are concerned, the relatively 
strong quantum reduction of the magnetization and the 
interplay between J t ss J s preclude a simple interpreta- 
tion of those two gapped phases based on the more subtle 
udd dimerized and uud trimerized ground states of the 
Ising-Heisenberg model. 

Last, the magnetization process of the spin-i Ising- 
Heisenberg and Heisenberg model on the second TIT 
lattice (Fig. [T]d) is displayed in Fig. |6] It is quite clear 
that the ground states of the other two Ising-Heisenberg 
and Heisenberg models are quite analogous with the ones 
previously described for the first TIT lattice (Fig. Q]i) 
and hence, there is no need to repeat the comprehen- 
sive discussion concerning the nature of possible ground 
states. It will be sufficient to mention the most signifi- 
cant differences. First, it should be pointed out that the 
ground states which have exactly the same character will 
eventually lead to the magnetization plateaus at different 
fractions of the saturation magnetization due to the fact 
that the number of quantum Heisenberg spins from the 
inner equilateral triangles is twice as large (the elemen- 
tary unit cell now consists of six spins from two inner 
equilateral triangles and one spin from the outer isosce- 
les triangle). Consequently, the spin-^ Ising-Heisenberg 
and Heisenberg models on the second TIT lattice ex- 
hibit the most robust magnetization plateaus at ^y, |, 
| and | of the saturation magnetization, which corre- 
spond to the udd dimerized state (fT2|) . the uud trimer- 
ized state (fT3|) , the uuu dimerized state (fT4| and the fer- 
rimagnetic state (|10|) . respectively. It is worth mention- 
ing that the classical (or quantum) ferrimagnetic state 
gives rise to a fractional magnetization plateau differ- 
ent from the uuu dimerized state (|14l) . which appear at 
| and | of the saturation magnetization, respectively. 
Owing to this fact, the spin-i Ising-Heisenberg model 
on the second TIT lattice shows in a small parameter 
range J t /J s G (±§ + f ) [or a G (70.1°, 73.3°)] a 

rather spectacular magnetization curve with four differ- 
ent intermediate plateaus, which reflect four consecutive 
field-induced transitions between the udd dimerized state 
(III), the uud trimerized state (IV), the uuu dimerized 
state (V), the classical ferrimagnetic state (I) and the 
saturated paramagnetic state (II) in a respective order 
along ascending field. 

To get a global view of the difference between the mag- 




FIG. 6: (a) A three-dimensional plot for the zero-temperature 
magnetization of the spin-i Ising-Heisenberg model on the 
second TIT lattice (Fig.[Tb) as a function of the magnetic field 
H/J and the parameter a; (b) The same three-dimensional 
plot but for the pure quantum spin-i Heisenberg model on 
the second TIT lattice of 21 sites. Note that since the number 
of spins is odd, the magnetization starts at f/2. 

netic behavior of the spin-i Ising-Heisenberg model on 
two considered TIT lattices, the relevant ground-state 
phase diagrams in the a — hj J plane are displayed in 
Fig. [7] for comparison. It is quite apparent from this 
figure that the classical ferrimagnetic state (I) generally 
extends over a much wider parameter region at the ex- 
pense of the udd dimerized (III), the uud trimerized (IV), 
and the uuu dimerized (V) ground states for the Ising- 
Heisenberg model on the second TIT lattice. In addition, 
it can be also clearly seen from Fig. [TJd that the magneti- 
zation curve of the Ising-Heisenberg model on the second 
TIT lattice indeed exhibits four consecutive field-induced 
transitions between the udd dimerized state (If 2[) . the uud 
trimerized state (JT3J), the uuu dimerized state (O, the 
classical ferrimagnetic state (| 10[) and the saturated para- 
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FIG. 7: The ground-state phase diagram of the spin-^ Ising- 
Heisenberg model on the first (Fig. [7^) and the second 
(Fig. [TJd) TIT lattice. The numbers in round brackets re- 
fer to the total magnetization normalized with respect to its 
saturation value. 

magnetic state (fTTj) . which are reflected in the relevant 
magnetization jumps between four different intermediate 
plateaus at Jj, ^, | and | of the saturation magnetiza- 
tion if a € (70.1°, 73.3°). A magnetization process of this 
type is not present in the Ising-Heisenberg model on the 
first TIT lattice, which may show at most three differ- 
ent intermediate plateaus before reaching the saturation 
magnetization (cf. Fig. [7^ with Fig. [7p). 



been studied with the help of numerical diagonalizations 
of finite-size clusters using the Lanczos algorithm. It has 
been demonstrated that the magnetization curves of the 
Ising-Heisenberg and Heisenberg models are quite remi- 
niscent even though the exact ground states for the Ising- 
Heisenberg model do not represent true ground states for 
the analogous quantum Heisenberg model. 

Our exact analytical results for the Ising-Heisenberg 
model shed light on an extraordinary diversity of the 
zero-temperature magnetization curves, which may in- 
clude up to three or four intermediate magnetization 
plateaus before reaching the saturation magnetization. 
As a matter of fact, we provide convincing evidence of 
the existence of three unconventional quantum ground 
states referred to as the udd dimerized state, the uud 
trimerized state and the macroscopically degenerate uuu 
dimerized state in addition to the classical ferrimag- 
netic and saturated paramagnetic ground states. The 
most robust magnetization plateaus of the quantum 
Heisenberg model remarkably coincide with the relevant 
magnetization plateaus of the Ising-Heisenberg model, 
which brings some insight into more complex quantum 
ground states of the Heisenberg model. Until now, 
the magnetic behaviors of the classical-quantum Ising- 
Heisenberg model and of the purely quantum Heisenberg 
model had only been compared for the particular case 
of the one-dimensional orthogonal-dimer chain^ The 
present study clearly demonstrates that a lot of insight 
into unusual quantum ground states can be achieved by 
looking at simpler classical-quantum spin models. 

Even though our analysis was restricted only to the 
zero-temperature magnetization process under which 
the most pronounced manifestations of quantum effects 
should be anticipated, we also expect interesting effects 
at finite temperatures. Indeed, a discrete lattice symme- 
try is broken in the ground state of all plateau phases, 
and we thus expect the system to undergo a phase transi- 
tion at finite temperature whenever the ground state is in 
a plateau state, as observed for instance in SrCu2 (603)2 
(see e.g. Ref . KLOh . The determination of the universality 
class of the transition for a given plateau, which should 
a priori depend on how the lattice symmetry is broken 
in this plateau, is a challenging problem that is left for 
future investigation. 



IV. CONCLUDING REMARKS 



In the present article, the ground state and zero- 
temperature magnetization process of the spin-^ Ising- 
Heisenberg and Heisenberg model on two geometrically 
frustrated TIT lattices have been investigated in de- 
tail. The magnetization process of the former classical- 
quantum Ising-Heisenberg model has been rigorously 
found using the eigenstates of the smallest commuting 
spin clusters. On the other hand, the magnetization pro- 
cess of the latter purely quantum Heisenberg model has 
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